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ABSTRACT 

Large-scale bulk flows in the Universe distort the initial density field, broadening the baryon- 
acoustic-oscillation (BAO) feature that was imprinted when baryons were strongly coupled to photons. 
Additionally, there is a small shift inward in the peak of the conventional overdensity correlation 
function, a mass-weighted statistic. This shift occurs when high density peaks move toward each 
other. We explore whether this shift can be removed by applying to the density field a transform 
(such as a logarithm) that gives fairer statistical weight to fluctuations in underdense regions. Using 
configuration-space perturbation theory in the Zel'dovich approximation, we find that the log-density 
correlation function shows a much smaller inward shift in the position of the BAO peak at low redshift 
than is seen in the overdensity correlation function. We also show that if the initial, Lagrangian den- 
sity of matter parcels could be estimated at their Eulcrian positions, giving a displaccd-initial-density 
field, its peak shift would be even smaller. In fact, a transformed field that accentuates underden- 
sities, such as the reciprocal of the density, pushes the peak the other way, outward. In our model, 
these shifts in the peak position can be attributed to shift terms, involving the derivative of the linear 
correlation function, that entirely vanish in this displaced-initial-density field. 
Subject headings: cosmology: theory - large-scale structure of Universe 



1. INTRODUCTION 



The "standard ruler" that baryon acoustic oscillations 
(BAO) provide is a major tool in cosmology. The BAO 
peak h as now been detected at unquestionable signifi- 
cance jElsenstein et al.||2005| |Cole et al.1|2005| |Ander-| 
son et al.|2012| ). The power arid cleanliness of the BAU 



method is such that despite these impressively large data 
sets, statistical uncertainties still dominate the system- 
atic effects. With the advent of even larger upcoming 
surveys, however, systematics such as the small shift in 
the BAO peak will soon be crucial to consider. 

The peak is distorted and shifted at low redshift, 
even in the conceptually straightforward real-space dark- 
matter field, without galaxy bias and redshift-space dis- 
tortions. This happens because of bulk displacements, 
understood most easily in a Lagrangian description. The 
broadening is largely removable by estimating the dis- 
placeme nt field with the simple, linear Zel'dovich approx- 



imation (lE iscnstcin et aL 2007 Padmanabhan et al. 2012 
Mehta et al. 2011 T Refinements in estimating the dis 



placement field, such as using the log density, are also 

|Lavaux et al.| [2 008 ; 
o|2U12| | Tassev &; Zal-j 



promising (Mohayaee et al. 



Kit 



2006 



aura 



& An 



•-•II 



Falck et al. 20Ti 

darriaga||2012[ ). 'While restoring the initial sharpness of 
the BAO peak has obvious statistical benefits, such re- 
construction methods can be rather involved. 

As surveys grow, the shift in the peak becomes more 
important to consider in analysis than the broadening. 
Various studies have attributed the shift to so-called 



"mode-coupling" terms in the power spectrum (Smith 


et al. 2008; Matsubara 2008; Sco et al. 


2008| |Padmanab- 


han & White||2009 Taruya et al.||2009 
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We can also understand the peak shift in configura- 
tion space. The overdensity correlation function £s( r ) — 
(5(x)5(r + x)) x is a volume average of the overdensity 
6 = p/p—1. The presence of the BAO peak in the linear 
correlation function means that on average in the ini- 
tial conditions, over- and underdensities are surrounded 
by over- and underdense shells, respectively, at a radius 
r = r s , the sound horizon. In the initial conditions, 
the Gaussian probability density function (PDF) of 8u n 
exactly balances the contribution of the over-dense and 
under-dense regions to £5. 

At low redshift, the PDF of 8 becomes highly skewed, 
as overdensities collect matter and become sharp peaks 
(where 8 grows arbitrarily high) at the expense of broad, 
shallow void regions (where 8 > — 1). If the BAO shell 
around each particle remained intact, tethered to move 
precisely with the particle toward or into clusters, the 
peak in £5 would not shift from its original location. 

However, 3D motions generally cannot preserve all 
BAO shells. Central over- and underdensities attract 
or repel their over- and underdense shells, broadening 
the averaged peak at low redshift. Moreover, the over- 
dense shells grow nonlinearly in density, and narrow 
(they must, from mass conservation); underdense shells 
broaden. This nonlinear density enhancement of over- 
dense shells at 8(r + x) gives more weight to <5(x) > 
than < 0, pulling the peak in £(r) to slightly lower r. 
This view of the pea k shift is similar to that of |Sherwin| 



& Zaldarriaga (2012 1, in which the peak shifts because of 
an increased abundance of tracers in overdense regions, 
where the universe expands more slowly than average. 

In this Letter, we use the Zel'dovich approximation 
(Zel'dovich 1970) to test the hypothesis that a den- 
sity transform that produces a more-Gaussian final- 
conditions density PDF leads to a reduction in the peak 
motion, as such a transform more equally weights the 
contributions in the volume average from overdense and 



2 



underdense regions. While we do not show mathemati- 
cally how transformations change the weightings in the 
correlation function average, which would involve ma- 
nipulating the PDFs of the transformed fields, we verify 
that the more weight low densities receive in a transfor- 
mation, the larger the peak position. We investigate two 
transformations that produce near-Gaussian PDFs - the 
log density and the displaced-initial-density - and one 
that over-weights low densities, (1 + S)^ 1 - 

Our investigation of log-density statistics is also moti- 
vated by recent work on clustering statistics of the log 
density and Gaussianized density, which have dramati- 
cally lower covariance on small sc ales, as well as bette r 
fidelity to the linear-theory shape ( Neyrinck et al.|2009 ). 
In principle, they can then gi ve much tighter constraints 



on cosmological parameters (Neyrinck 20111. In fact, 
analysis of the log-density correlation iimction accesses 
all of the Fisher information in an idealized entirely log- 
normal field, whereas using even arbitrarily high-point 
S correlation functions can give only a small fraction of 
that total information ( Carron|2011 Carron & Neyrinck 
2012 1. Here, for the first time we explicitly show the ben- 
efits of a log transform on large scales, where the benefits 
might have been thought to be negligible. 

2. THEORY 

In this section, we derive the nonlinear correlation 
functions of the Zel'dovich density, the log density, and 
the displaced-initi al-density field. We follow the ap- 
proach outlined by McCullagh & Szalay (2012) for per- 
turbatively calculating the nonlinear density correlation 
function in the Zel'dovich approximation. We briefly re- 
view this approach below, and extend it to calculate the 
nonlinear correlation functions of the log density and the 
Lagrangian density 5^ at an Eulerian position x. For the 
latter, we examine the effect of Zel'dovich displacements 
by themselves, without considering the change in density 
at a constant Lagrangian coordinate. 

We start with the Zel'dovich approximation, which 
maps particles' initial Lagrangian coordinate, q, to their 
co-moving Eulerian coordinate, x, through the gradi- 
ent of the initial displ acement potent ial, 4>(q), and the 
growth function D(t) ( |Zel'dovich||l970[ ). The initial dis- 
placement potential is related to the initial density field 
through the Poisson equation: 



x(q,t) = q-D(t)V q (f>(q) 
D(t)V^(q) = <5 Q (q,*) . 



(1) 



The density can be written in terms of the Jacobian of the 
transformation from Lagrangian to Eulerian coordinates 
from the conservation of mass. 



p(x.,t) 



dxj 
dq 3 



1 



J(q,t) 



1 + *(q(x)) . (2) 



Here, 8 is the nonlinear overdensity as a function of La- 
grangian position. Equation (pi) for the Eulerian density 
is only strictly valid before shell-crossing, where the map- 
ping from q to x is one-to-one. 

The Jacobian can be written in terms of the eigen- 
yalues (Ai, A 2 , a nd A3) of the symmetric matrix dy(q) 
( Zerdovich||1970 ), or equivalently in terms of the invari- 
ants of the matrix (ii, I2, and I3), which are simple 



functions of the eigenvalues ( Bouchet et al.|[T995 ) 
d 2 Mq) 



J(q,t) 



dqidqj 

(1-DX 1 )(1-DX 2 )(1-D\ 3 ) 
1 - D/i(q) + D 2 I 2 {q) - D 3 I 3 {q) 



(3) 

(4) 
(5) 



The nonlinear overdensity can thus be expressed in 
terms of the initial quantities by a Taylor expansion of 
the inverse Jacobian, to any order: 

S(q,t) = Dh(q) + D 2 (^(q) 2 - / 2 (q)) 

+ D 3 (^(q) 3 - 27 1 (q)7 2 (q) + J 3 (q)) + ... 

(6) 

In order to express the overdensity as a function of 
final (Eulerian) position, x, we Taylor expand Equation 

§ about the point x = q. The nonlinear density as a 
iction of x to third order in D is then: 

d 2 4>{q) d4>(q) 8S(q,t) 



dqidqj dqj dqi 



l D 2y d 2 S(g 7 t) d^g) 90(g) 
2 4-^ dqidqj dqi dqj 



q=x 



(7) 

The correlation function in co-moving Eulerian coor- 
dinates can be written in powers of D using the above 
expansion of the Eulerian overdensity in terms of the 
initial quantities. Because the initial overdensity is as- 
sumed to be a zero-mean Gaussian random field, the odd 
moments vanish. The first two terms of the correlation 
function are then: 

(S (x, t)*(x + r, t)> = ^ (r)D 2 + (r)# 4 



&(r,t) 



We define the functions: 
1 



C(0 = 



2tt 2 



PUk)] n (kr)k m+z dk 



(8) 



(9) 







where j n is the spherical Bessel function of order n and 
Pn n (fc) is the power spectrum of the initial density fluc- 
tuations. Using this definition, the linear term in the 
expansion of the correlation function is: 

&Hr)=e Q (r), (10) 

the spherically symmetric Fourier transform of the linear 
power spectrum. 

The first nonlinear term in the correlation function can 
be calculated using sph erical harmonics. For more d etails 
on the calculation, see McCullagh & Szalay (2012). The 
final expression is: 

£f (r) = ^°(rf + ^(rf + yi{rf - ^(rj&r) 



(11) 
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We can use the same approach to calculate the cor- 
relation function of the log density field. Starting with 



Equation 
Equation 
define as t 



We label the correlation function of this quantity as 
£s L ■ The first two terms are: 



2|) and using the definition of the Jacobian in 
we write the log of the density, which we 
le quantity -A(q), as: 



In (1 + 6{q, t)) = A(q) = - In J(q, t) 

= -{\n{l - DXi) + ln(l - D\ 2 ) 

+ ln(l-£>A 3 )) . (12) 

We expand this expression in powers of D, and rewrite 
it in terms of the invariants of the deformation tensor: 



d>) = £oV) 

^ ) M = -k(0Ko o M + k o M : 



i 



Co~ 8 (0)^(r) + tj£ 



2 

+ 3 

-2/ \ e 2 



C 2 °(r) 2 
(r)Cf(r) 



1 



(19) 
(20) 



Again, we note the damping term, — 3C0 (0)Co ( r ): 
which is the same as we found in £4. We also note that 



ln(l + <f(q,t)) 



D(A a + A 2 + A 3 ) + l -D 2 {\\ + A 2 + A3) 



the expression for ^ 2 -* is simpler than either ^ 2 -* 



t(2) 

or^V- 



+ 3^ 3 (A? + A3 + A 3 ; ) + 



I>Jl 



2 1 1 
3 v 1 



21, 



3hh + 3/3) 



(13) 



(14) 



Finally, to investigate the peak shift in a transforma- 
tion which further boosts the weight of underdensitics, 
we consider the correlation function of the reciprocal of 
the density, £i/ p - This statistic is of some theoretical in- 
terest in a Lagrangian approach, since (1 + (5) _1 is simply 
the Jacobian, Equation Q. When we expand to Eule- 
rian coordinates we get: 



We now have the equivalent of Equation but for 
the log density field. We transform this into Eulerian 
coordinates as above to get an expression for A(x,t). 
Because this quantity has a non-zero mean, we define 
the correlation function as: 



(l + aM)-^ (V(q) + £ 2 ]r 



^ 3 E 



dcj>(q) dJ(q) 
dqi dq t 



d 2 0(q) afr(q) 0J(q) 
dqidqj dqj dqt 



U(r, t) = ((A(x, t) - A)(A( X + r, t) - A)) 



(15) 



where A is the Eulerian mean. 
The first two terms in the correlation function are then: 



D 3 y d 2 J(q) d<j)(q) d^q) 
2 ^-i dqidq dq t dq 3 



(16) 



The first two terms of £i/ p are: 



13 



20 



r(!) 



30 



C, '(r-KJ(r) 

^ 2 (o)e 2 w 



2 
3 



21 

_1 (r)^(r) 



35 



q=x 

(21) 



(22) 



(17) 



Note that the first term in £ A (r) can be written as 

— |o"o^A ( r )> wriere o 2 is the variance of the initial den- 
sity field (where the density is assumed to be smoothed 
on some scale). This term describes the reduction in 
the amplitude of £,a compared with the linear correla- 
tion function. 

Next, we look at the initial density field at Eulerian 
position x. For this quantity, we take the initial density, 
(5o(q), and use the Zel'dovich formula, Equation (JlJ, to 
relate q to x. The expression for <$l(x) (where the sub- 
script L stands for Lagrangian) is equivalent to Equation 
Q, but with S replaced with <5 : 



§6TV)£(r) - ^3»? 3 V) + ^o~V)£o» 



(23) 



<5 L (x) = 6 (q)+D 



2 E 



D 



3 E 



o 2 



d^(q) dS (q) 
dqi dqi 

q) 00(q) 8S (q) 



2 



E 



dqidqj dq 3 dq t 
d 2 6 (q) dcb(q) d<j>(q) 



Because our interest in this transformation is to test 
a more extreme weighting of underdensities, and we do 
not anticipate a practical measurement of £i/ p (which, 
for example, could be violently sensitive to discreteness 
noise), for clarity we do not plot £i/ p in the figures below, 
but do discuss its peak position. 

3. NUMERICAL RESULTS 

Our goal is to understand how the various transforma- 
tions to the density field affect the correlation function, 
specifically at the BAO peak. In this section, we con- 
sider the effects of the first nonlinear term in the corre- 
lation functions on the BAO peak position and shape. 
We then examine the various terms in the expressions 
in order to understand their role in the overall correla- 
tion function. For the following we use an initial p ower 



spectrum generated by CAMB flLewis et al. 20001, as 



dqidqj dqi 



(18) 



suming fiducial WM AP5 cosmological parameters (Ko- 



( |Lewis 
)logica 



matsu et al. 2009): [H Q , fi A , a 



08 ) 



CDMj ^b, n, 

(70.1, 0.233, 0'.tHoT 0.721, 0.96, 0.817). We smooth the 
linear power spectrum with a a — 5 Mpc//i Gaussian 



4 



which was a large enough smoothing length to give agree- 
ment between the analy tical expression and numer ical 
Zel'dovich realizations in McCullagh & Szalay (2012) up 
to z = 0. 

Fig. [T] shows the correlation functions at z — for 
the linear theory density, Zel'dovich density, Zel'dovich 
log-density, and displaced-initial-density. We note that 
the BAO peak of all nonlinear correlation functions are 
broadened compared with £i; n . The smoothing of the 
acoustic peak relative to that of the linear correlation 
function is due to large-scale bulk motions, which are 
present in the three nonlinear fields. In the case of £g we 
also observe an enhancement at the peak, perhaps from 
over-weighting overdense mass elements, which contract. 
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Fig. 1. — Nonlinear correlation functions at z = 0. Green dashed 
line - linear theory, Blue dotted line - nonlinear density from 
Zel'dovich approximation, magenta solid line - displaced-initial- 
density field, 5^ (x), black dot-dashed line - ln(l+<5) from Zel'dovich 
approximation 



Fig. [2] shows a close-up of the BAO peak, and the po- 
sition of each of the peaks (r p ) is indicated. We mea- 
sure the peak of r 2 £ because this is the quantity that 
is integrated to give the variance on a given scale. The 
peaks in order of largest to smallest r p are: £i/ p : r p = 
106.5 Mpc/h (not shown), £ lin : r p = 106.4 Mpc/h, 



6 



= 106.2 Mpc/h, £4 



106.1 Mpc/h, and 



= 105.8 Mpc/h. There are other roughly equiva- 



lent definitions of the peak position, such as the minimum 
X 2 of template correlation functions or the peak through 
a wavelet filter. These are essential to use in the pres- 
ence of cosmic variance, but the simpler peak definition 
used here is acceptable since we use a linear correlation 
function with no cosmic variance. 

Next, we look at the effects of individual terms in the 
nonlinear expressions on the location of the peak. Fig. 
[3] shows the terms that contribute to each of these non- 
linear correlation functions. We exclude the damping 
terms of £4 and t;$ L because they do not contribute to 
the shift in the peak position. Each term is represented 
by a different color, and the solid lines have a different 
amplitude in each expression. The dot-dashed line has 
the same amplitude in £5 and £4 (and is zero in £g L ) , and 
the dashed line is the same in all three expressions. From 
this plot we see that the effect of the dashed black line 
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Fig. 2. — Peak locations of the various correlation functions 
at 2 = 0. Green dashed line - linear theory, peak is at r = 
106.4Mpc//i. Blue dotted line - nonlinear density from Zel'dovich 
approximation, peak is at r = 105.8 Mpc/h. Magenta solid line - 
displaced-initial-density, &l(x), peak is at r = 106.2 Mpc/h. Black 
dot-dashed line - ln(l + 8) from Zel'dovich approximation, peak is 
at r = 106.1 Mpc/h. 

is to smooth the peak in all three cases, by increasing 
the amplitude around the peak and decreasing the am- 
plitude at the peak. Several of the other terms cause a 
shift in the peak position by showing a slope around the 
peak. For £$ L , the magenta line is the only term with this 
property, so the shift is not very large. In £,$ and £4, the 
red (dot-dashed) and magenta (solid) lines have nearly 
opposite slopes, and so the shift from the combination of 
these terms is small. 

The major contributor to the shift in the peak position 
in £4 and £5 is the (solid) blue line, which is the term 

proportional to — ( r )Cf 1 ( r )- I n £<5 the amplitude of this 
term is 16/5, whereas for £4 the amplitude is 6/5. Thus 
we sec a larger shift in the peak position in £5 than in 
£a> and both have larger shifts than £g L . We also note 
that in the expression for £i/ p this term has the opposite 
sign, indicating a shift to greater radius. 

We can understand this term better by using the re- 
cursion relation for spherical Bessel functions: 



ioO) = ~h{x) 



(24) 



This allows us to rewrite £™(r) in terms of derivatives of 
the linear correlation functions with respect to r: 



^KrV) = £o°MVM' 



(25) 
(26) 



In fac t this is exactly the term [Crocce &: Scoccimarro| 
(2008) found to be responsible for the BAO peak shift 
from the mode coupling term of the correlation function 
in renormalized perturbation theory. This term is the 
product of the derivative of the linear density correlation 
function (£o( r )) an d the derivative of the linear velocity 
correlation function {£,Q 2 {rj). 

4. CONCLUSION 

In this Letter, we have explored the idea that the shift 
in the position of the BAO peak in the usual density cor- 
relation function is due to over- weighting of high density 
peaks that bulk flows bring toward each other. We have 
shown that the two transformations on the density field 
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Fig. 3. — Top to bottom: Terms in the expressions for £ 



(2) ("21 

£ 4 , and Q J (not including damping terms in the last two). Dif- 
ferent terms are represented by different colors: Black dashed - 
terms with same amplitude in all 3 expressions (l/3(£Q(r)§^~ 2 (r) — 
?o( r )?r7 2 ( r )) + 2 / 3 €l( r )^2 2 ( r ))' Red dot-dash - terms that are the 
same in both & and £ A (-4/5^(r)^ _1 (r) + 4/35^(r) 2 ). The am- 
plitude of this term is zero for £g L . Green - amplitude of £o( r ) 2 



term. Magenta - amplitude of ^(r) 2 term. Blue - amplitude of 

eiO-^rV). 

that suppress overdense regions and boost underdense re- 
gions in statistical weight result in a reduced shift in the 
BAO peak position at low redshift. Moreover, we have 
shown that a transformation that further weights under- 
dense regions causes the peak to shift in the opposite 
direction, to larger r. We used the Zel'dovich approx- 
imation and configuration-space perturbation theory to 
examine the first nonlinear contribution to the correla- 
tion functions of the nonlinear density, the log-density, 



and the displaced-initial-density fields. We identified the 
terms largely responsible for the shift in the BAO peak, 
and found that they were proportional to the first deriva- 
tive of the linear correlation function. The displaced- 
initial-density exhibits the smallest shift due to a lack of 
shift terms in the nonlinear correlation function expres- 
sion. 

The log-density correlation function has some shift, 
but it is smaller than in the nonlinear density corre- 
lation function. In fact, the correlation function of 
the log density is almost indistinguishable from that of 
the displaced-initial-density on BAO scales. The near- 
equivalence of the log density and the displaced initial 
densit y recalls the lognormal model of IColes & Jones 
( |1991| , in which gravity enhances initia 



densities with 



an exponential transform, even though a Zel'dovich- 
realization density field is not an exactly lognormal field. 

By studying both the nonlinear density and displaced- 
initial-density, we decouple the effects of density en- 
hancements from the bulk flows of initial fluctuations. 
The fundamental broadening of the initial peak is seen 
in £g L , whereas in {;$ the peak is additionally sharpened. 
However, we suspect that this sharpness is not statisti- 
cs) cally useful because it comes from only the fraction of 
the mass elements that have contracted. We intend to 
test the noise properties of the transformed fields in N- 
body simulations. Future work is also needed to exam- 
ine the potential of transformations such as these in the 
presence of galaxy bias, discreteness, and redshift-space 
distortions. 
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